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Abstract 

We explore the effects of social influence in a simple market model in which a large number 
of agents face a binary choice: 'to buy/not to buy' a single unit of a product at a price posted 
by a single seller (monopoly market). We consider the case of positive externalities: an agent 
is more willing to buy if the other agents with whom he/she interacts make the same decision. 

We compare two special cases known in the economics literature as the Thurstone and 
the McFadden approaches. We show that they correspond to modeling the heterogenity in 
individual decision rules with, respectively, annealed and quenched disorder. More precisely 
the first case leads to a standard Ising model at finite temperature in a uniform external field, 
and the second case to a random field Ising model (RFIM) at zero temperature. 

We illustrate some dynamic properties of these models, making use of numerical simula- 
tions in an ACE (Agent based Computational Economics) approach, and we study analytically 
the equilibrium properties in the mean field limit. 

Considering the optimisation of profit by the seller within the McFadden/RFIM model 
in the mean field limit, we exhibit a new first order phase transition: if the social influence 
is strong enough, there is a regime where, if the mean willingness to pay increases, or if the 
production costs decrease, the optimal solution for the seller jumps from one with a high price 
and a small number of buyers, to another one with a low price and a large number of buyers. 

This paper has been presented at the interdisciplinary conference WEHIA 2003: it is 
intended to be accessible to economists and physicists - without, hopefully, discouraging any 
of them! 

1 Introduction 



Following Kirman jlUl lllj , we view a market as a complex interactive system with a 
communication network. We consider a market with discrete choices 0, and explore 
the effects of localised externalities (social influence) on its properties. We focus on 
the simplest case: a single homogeneous product and a single seller (monopoly). On 
the demand side, customers are assumed to be myopic and non strategic. The only 
cognitive agent in the process is the monopolist, who determines the price in order 



to optimise his profit. We obtain interesting and complex phenomena, that arise 
due to the interactions between heterogeneous agents. 

Agents are assumed to have idiosyncratic wilhngness-to-pay (IWP) which are 
described by means of random variables. We consider two different models: on one 
hand, the IWP are randomly chosen and remain fixed, on the other hand, the IWP 
present independent temporal fluctuations around a fixed (homogeneous) value. We 
will show that they correspond to statistical physics models, in the former case with 
quenched disorder, and in the latter with annealed disorder. In both cases we assume 
that the heterogeneous preferences of the agents are drawn from the same (logistic) 
distribution. The equilibrium states of the two models generally differ, except in 
the special case of homogeneous interactions with complete connectivity. In this 
special situation, which corresponds to the mean-field approximation in physics, the 
expected aggregate steady-state is the same in both models. 

In the following, we first present the demand side, and then consider the op- 
timisation problem left to the monopolist, who is assumed to know the demand 
model and the distribution of the IWP over the population, but cannot observe the 
individual (private) values. 

2 Simple models of discrete choice with social influence 

We consider a set iljsr N agents with a classical linear IWP function Each 
agent i £ either buys (wj = 1) or does not buy {uji = 0) one unit of the single 
given good of the market. A rational agent chooses uJi in order to maximise his 
surplus function Vi : 

max Vi= max uJi{Hi + JikOJk - P), (1) 
<.^,e{o,i} ^,G{o,i} 

where P is the price of one unit and Hi represents the idiosyncratic preference 
component. Some other agents k within a subset "di G (such that k € t?i) 
hereafter called neighbours of i, influence agent i's preferences through their own 
choices cok- This social influence is represented here by a weighted sum of these 
choices. Let us denote Jik the corresponding weight i.e. the marginal social influence 
of the decision of agent /c G i?j on agent i. When this social influence is assumed to 
be positive (Jj^ > 0), it is possible, following Durlauf (II, to identify this external 
effect as a strategic complementarity in agents' choices 3_. 

For simplicity we consider here only the case of homogeneous influences, that 
is, identical positive weights Jik = and identical neighbourhood structures of 
cardinal n, for all the agents. That is, 

Jik = Ji) = J/n > yiG^N, k £ 'di (2) 

2.1 Psychological versus economic points of view 

Depending on the nature of the idiosyncratic term Hi, the discrete choice model 
may represent two quite different situations. Following the typology proposed by 
Anderson et al. |Tj, we distinguish a "psychological" and an "economic" approach 
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to individual choices. Within the psychological perspective (Thurstone |2H|)) the 
utility has a stochastic aspect because "there are some qualitative fluctuations from 
one occasion to the next... for a given stimulus" (this point of view will be referred 
to hereafter as the TP-model). On the contrary, for McFadden U2] each agent has 
a willingness to pay that is invariable in time, but may differ from one agent to 
the other. In a "risky" situation the seller cannot observe each specific IWP, but 
knows its statistical distribution over the population (we call this perspective the 
McF-model). Accordingly, the TP and the McF perspectives differ in the nature of 
the individual willingness to pay. 

In the TP model, the idiosyncratic preference has two sub-components: a con- 
stant deterministic term H (the same for all the agents), and a time- and agent- 
dependent additive term ei{t), so that Hi = H + €{. The ei{t) are i.i.d. random 
variables of zero mean; in the simulations they are refreshed at each time step 
(asynchronous updating). Agent i decides to buy according to the conditional prob- 
ability 

P{uJi = l\z,{P,H))=V{e, > Zi{P,H)) = l-F{zi{P,H)), (3) 

with 

Zi{P,H) = P-H~J^ Y^uJk, (4) 

where F{zi) = V{ei < Zi) is the cumulative distribution of the random variables e^. 
In the standard TP model, the agents make repeated choices, and the time varying 
components ei{t) are drawn at each time t from a logistic distribution with zero 
mean, and variance = 7r^/(3/3^): 

^<-' = l + exp(-/?.) - 

In the McF model, the private idiosyncratic terms Hi are randomly distributed 
over the agents, but remain fixed during the period under consideration. There are 
no temporal variations: the are strictly zero. In analogy with the TP model, it 
is useful to introduce the following notation: Hi = H + 9i, and to assume that the 
6i are logistically distributed with zero mean and variance = 7r^/(3/3^) over the 
population. This assumption implies: 

hm 4 y = and lim ^YHi = H. (6) 

Thus, the correspondence between models is better the larger the number of agents, 
and is strict only in the limit of an infinite population. In practice, the simulations 
presented in the following sections show that the theoretical predictions are already 
verified for population sizes of the order of some hundreds. 

For a given distribution of choices in the neighbourhood 'di, and for a given price, 
the customer's behaviour in the McF model is deterministic. An agent buys if: 

ei>P-H-J^ Y^ojk. (7) 
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2.2 'Annealed' versus 'quenched' disorder 



Since we have assumed isotropic (hence symmetric) interactions, there is a strong 
relation between these models and Ising type models in Statistical Mechanics, which 
is made explicit if we change the variables coi G {0, 1} into spin variables Si G {±1} 
through uji = All the expressions in the present paper can be put in terms 

of either Sj or LOi using this transformation. In the following we keep the encoding 
uJi £ {0,1}. The assumption of strategic complementarity in economics corresponds 
to having ferromagnetic couplings in physics (that is, the interaction J between Ising 
spins is positive). In this case, the spins Si (and consequently the choices cjj) all tend 
to take the same value. 

In physics, the TP model corresponds to a case of annealed disorder. Having a 
time varying random idosyncratic component is equivalent to introducing a stochas- 
tic dynamics for the Ising spins. In the particular case where F(z) in © is the 
logistic distribution, we obtain an Ising model in a uniform (non random) external 
field H — P, at temperature T = The McF model has fixed heterogeneity; it 
is analogous to a Random Field Ising Model (RFIM) at zero temperature, that is, 
with deterministic dynamics. The RFIM belongs to the class of quenched disorder 
models: the values Hi are equivalent to random time-independent local fields. The 
"agreement" among agents favored by the ferromagnetic couplings may be broken 
by the influence of these heterogeneous external fields Hi. Due to the random distri- 
bution of Hi over the network of agents, the resulting organisation may be complex. 
Thus, from the physicist's point of view, the TP and the McF models are quite 
different: uniform field and finite temperature in the former, random field and zero 
temperature in the latter. The properties of disordered systems have been and still 
are the subject of numerous studies in statistical physics. They show that annealed 
and quenched disorder can lead to very different behaviours. 

The TP model is well understood. In the case where the agents are situated on 
the vertices of a 2-dimensional square lattice, and have four neighbours each, there 
is an exact analysis of the model for P = P„ (the unbiased case, . see below) 
due to Onsager |[14j. Even if an analytical solution of the optimization problem 

for an arbitrary neighbourhood does not exist, the mean field analysis gives 
approximate results that become exact in the limiting situation where every agent 
is a neighbour of every other agent {i.e. all the agents are interconnected through 
weights (HJ). On the contrary, the properties of the McF model are not yet fully 
understood. A number of important results have been published in the physics 
literature since the first studies of the RFIM by Aharony and Galam [SI [Tj (see also 

j^lj)- Several variants of the RFIM have already been used in the context of 
socio-economic modelling jHl El 121 121] • 

2.3 Static versus dynamic points of view 

Hereafter, we restrict our investigation to the McF model in the case of a "global" 
externality. That is, we consider homogeneous interactions and full connectivity, 
which is equivalent to the mean field model in physics. Within this general frame- 
work, we are interested in two different perspectives. First we consider a static point 
of view: solving for the equality between demand and supply we determine the set 
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of possible economic equilibria. This will allow us to analyse in section^the optimal 
strategy of the monopolist, as a function of the model parameters. 

Next (section we consider the market's dynamics assuming myopic and non 
strategic agents: based on the observation of the behaviour of the other agents at 
time t — 1, each agent decides at time t to buy or not to buy. We show that, in 
general, the market converges towards the static equilibria of the preceding section, 
except for a precise range of the parameter values where interesting static as well as 
dynamic features are observed. 

These two kinds of analysis correspond in Physics to the study of the thermal 
equilibrium properties within the statistical ensemble framework on the one hand, 
and the out of equilibrium dynamics (which, in most cases, approaches the static 
equilibrium through a relaxation process) on the other hand, respectively. 



3 Aggregate demand 

As discussed in the preceding section, we consider the full connectivity case with 
isotropic interactions given by © with cardinal n = — 1, in the limit of a very 
large number of agents. In this limit the penetration rate r], defined as the fraction 
of customers that choose to buy at a given price: rj = limj\f^oo J2k=i ^k/^ ^ can be 
approximated by Jlke-d^k/i^ — 1), which is proportional to the social influence 
term of the agents' surplus function (eq. (^). Equation ^ may thus be replaced 
by 

ei> P - H - J r]. (8) 

For the following discussion it is convenient to identify the marginal customer, indif- 
ferent between buying and not buying. Let Hm = H + 9m^^ his IWP. This marginal 
customer has zero surplus {Vm = 0), that is: 

em = P-H-Jrj = z. (9) 

Thus, an agent i buys if 6i > 6m, and does not buy otherwise. Equations Q and 
allow us to obtain 77 as a fixed point: 

r] = l- F{z) (10) 

where z, defined in © depends on P, H, and rj, and F{z) is the cumulative distribu- 
tion of the IWP around the average value H. Note that this (macroscopic) equation 
is formally equivalent to the (microscopic) individual expectation that = 1 in the 
TP case ©. Using the logistic distribution for 6i, we have: 

_ 1 

l + exp(+/3z) 

Equation ()1U|) allows us to define r] as an implicit function of the price through 

$(r/, P) = 7]{P) + F{P -H-J 7]{P)) -1 = 0. (12) 

The shape of this (implicit) demand curve can be evaluated using the implicit deriva- 
tive theorem: 

dri{P) -d^/dP -f{z) 



dP d^/dr] l-Jf{z) 



(13) 



5 



where f{z) = dF{z)/dz is the probabihty density. 

Since for a given P, equation (jllj) defines the penetration rate r/ as a fixed-point, 
inversion of this equation gives an inverse demand function: 

pd{ri)=H + J (14) 
P V 

At given values of /3, J and H, for most values of P, (|11|) has a unique solution 
Ti{P). However for (3 J > (3Jb = 4, there is a range of prices 

Pi{(3J, (3H)<P< P2ipJ, (3H) (15) 

such that, for any P in this interval, has two stable solutions and an unstable 
one. The limiting values Pi and P2 are the particular price values obtained from 
the condition that eq. (|11|) has one degenerate solution: 

„ = l-f(.),and 

dr] 

This gives (3Jrj{l — ry) = 1 (together with eq. (fTT|) ). This equation has two solutions, 
m < 1/2 < r?i, 




; {1,2} (16) 



Note that r]i depends only on (3J. Then, the limiting prices Pj are equal to the 
inverse demands associated with these values T]i, which are, from (|14j) : 

P^ = H + Jr^, + \ H^—^] ; i G {1, 2}. (17) 

p m 

Note that these limiting prices are not necessarily positive. 

It is interesting to note that the set of equilibria is the same as what would be 
obtained if agents had rational expectations about the choices of the others: if every 
agent had knowledge of the distribution of the Hi , he could compute the equilibrium 
state compatible with the maximisation of his own surplus, taking into account that 
every agent does the same, and make his decision (to buy/not to buy) accordingly. 
For PJ < (3Jb every agent could thus anticipate the value of rj to be realized at 
the price P, and make his choice according to ©• For f3J > pJs, however, if the 
price is set within the interval [Pi,P2], the agents are unable to anticipate which 
equilibrium will be realized, even though the one with the largest value of rj should 
be preferred by every one (it is the Pareto dominant equilibrium). 



4 Supply side 

On the supply side, we consider a monopolist facing heterogeneous customers in a 
risky situation where the monopolist has perfect knowledge of the functional form 
of the agents' surplus functions and the maximisation behaviour He also knows 
the statistical (logistic) distribution of the idiosyncratic part of the reservation prices 
{Hi). In the special case of global externality considered here, where the interactions 
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are the same for all customers, as in equation (jSJ with n = N—1, as just seen, the TP 
model and the McF one have the same equilibrium states. Indeed, the monopolist 
cannot observe any individual reservation price. He observes only the cumulated 
result of the individual choices (which are either to buy or not to buy). As a result, 
the conditional probability for an agent taken at random by the monopolist to be 
a customer is formally equivalent in both cases. Accordingly, hereafter we limit 
ourselves to the case where the demand follows the McF model, in the limit of full 
connectivity. The social influence on each individual decision is then close to J ry, 
and the fraction of customers rj is observed by the monopolist. That is, for a given 
price, the expected number of buyers is given by equation (|Tn|) . 



4.1 Profit mcLximisation 

Let C be the monopolist cost for each unit sold, so that 

p = P-C (18) 
is his profit per unit. Since P — H = {P — C) — {H — C), defining 

h = H -C, (19) 

we can rewrite z in (j^ as: 

z = p-h-Jr]. (20) 

Hereafter we write all the equations in terms of p and h. 

Since each customer buys a single unit of the good, the monopolist's total ex- 
pected profit is p N 7]. Thus, in this mean field case, the monopolist's profit is 
proportional to the total number of customers. He is left with the following max- 
imisation problem: 

PM = argmaxn(p), (21) 

p 

where N IL{p) is the expected profit, with: 

U{p)=pr]ip), (22) 

and r]{p) is the solution to the implicit equation (|1()|) . If there is no discontinuity in 
the demand curve r]{p) (hence for pj < 4), pM satisfies dli{p)/dp = 0, which gives 
drj/dp = —rj/p at p = pM- Using the implicit derivative \\\\\. we obtain at p = pM- 

■^(^) - a (23) 

where z, defined in (^OJ, has to be taken at p = pM- 

Because the monopolist observes the demand level r], we can use equation (|lflj) 
to replace 1 — F{z) by r]. After some manipulations, equation (|^^ gives an inverse 
supply function p^{ri): 
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We obtain pM and rjM as the intersection between supply (|2Hl and demand (fTH) : 



PM=p'{ilM)=P'^{md), (25) 

where p'^ = P'^ — C. 

The (possibly local) maxima of the profit are the solutions of (|25|) for which 

^ < »■ P") 

It is straightforward to get the expression for the second derivative of the profit: 

2?7 - 1 



dp^ p 



1 + 



2/3p(l - r^f 



(27) 



from which it is clear that the solutions with r/ > 1/2 are local maxima. For rj < 1/2, 
condition H26(l reads 

1-2?? 

2/3p(l - ^)2 < ^^^^ 
Making use of the above equations, this can also be rewritten as 

2(3Jr]{l-r]f<l. (29) 

For PJ > 13 J B = 4, the monopolist has to find p = pm which realises the 
programme: 

PAf :max{n_(p^),n+(p*^)} (30) 

p^ = ar (7 max n_|_ (p) = p r]^{p), (31) 

p^ = arg max n_ {jp) = p rj- (p) (32) 

where the subscripts + and — refer to the solutions of (|11|) with a fraction of buyers 
larger, respectively smaller, than 1/2. 

To illustrate the behaviour of these equations, figure ^ represents several exam- 
ples of inverse supply and demand curves corresponding to different market config- 
urations. 



4.2 Phase transition in the monopohst's strategy 

In this section we analyse and discuss the solution of the optimal supply-demand 
static equilibria, that is, the solutions of equations and (jSHJ- As might be 
expected, the result for the product (3pM depends only on the two parameters (3h 
and /3J. That is, the variance of the idiosyncratic part of the reservation prices fixes 
the scale of the important parameters, and in particular that of the optimal price. 

Let us first discuss the case where h > 0. It is straightforward to check that in 
this case there is a single solution rjM- It is interesting to compare the value of pM 
with the value Pn corresponding to the neutral situation on the demand side. The 
latter corresponds to the unbiased situation where the willingness to pay is neutral 
on average: there are as many agents likely to buy as not to buy (ry = 1/2). Since 
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Fi gUrG 1: Inverse supply and demand curves p^{Tj) and p^{r})^ for different values of h and J (/3 = 1, 
C = hence h = H). The equilibrium prices are obtained at the intersection between the demand (black) 
and the supply (grey) curves. 

The two graphics on top illustrate the difference between a complete absence of externality (J = 0, dashed 

lines) and a strong externality (J = 4, solid lines). The case h = 2 (left) corresponds to a strong positive 

average of the population's IWP (h = 2), whereas the population is neutral for h = (right). 

The four graphics labelled (a) to (d) were obtained for values of h and J corresponding to the points (a) 

to (d) in the phase diagram (figure |2J. They all have negative values of the average of the population's 

IWP {h < 0), so that in the absence of externality only few consumers would be interested in the single 

commodity. 

(a) corresponds to the coexistence region between two local market equilibria in figure El but one of them 
(not shown) is not relevant since it corresponds to a negative price solution, (b) lies in the coexistence region; 
in this case, the optimal market equilibrium is the one with high rj. (c) lies in the region with only one 
market equilibrium, with few buyers (small rj). (d) corresponds to a large social effect; the single market 
equilibrium has large r] and a high price. 
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Figure 2; Fraction of buyers rj, optimal price f3pM and monopolist profit Hm, as a function of the 
social influence, for /3/i = —2. The superscripts — and + refer to the two solutions of equations I25i that 
are relative maxima. 




Figure 3: Phase diagram in the plane {(3 J, Ph}: the grey region represents the domain in the parameter 
space where coexist two maxima of the monopolist's profit, a global one (the optimal solution) and a local 
one. Inside this domain, as /3J and/or (3h increase, there is a (first order) transition where the monopolist's 
optimum jumps from a high price, low penetration rate solution = r;_ to one with low price, large r] = rjj^. 
The circles on the transition line have been obtained numerically, the smooth curves are obtained analytically 
(see the Appendix and 1131 for details). The points (a) to (d) correspond to the inverse supply and demand 
curves represented in figure ^ In the white region, for /3J < 27/8, the fraction of buyers, r], increases 
continuously from to 1 as /3h increases from — oo to +oo (c-d). At the singular point A, (jSJ = 27/8, 
0h = —3/4 — log(2)), r)-|- = -q— = 1/3. At point B (/3J = 4, (3h = —2), the local maximum with t) large 
appears with a null profit and rfj^ = 1/2. In the dark-grey region below B, this local maximum exists with 
a negative profit, being thus non viable for the monopolist (a). 
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the expected willingness to pay of any agent i \s h + 6i + J /2 — its average over 
the set of agents is /i + J/2 — p. Thus, the neutral state is obtained for 



Pn = h + J/2. (33) 

To compare pM with it is convenient to rewrite equation 1)14(1 as 

Pip'' - Pn) = PJ(.r] - 1/2) + ln[(l - 7?)/7?]. (34) 

This equation gives p'^ = Pn for rj = 0.5, as it should. For this value of rj, equation 
(PH) gives p^ = Pn only if + J) = 2: for these values of J and h, the monopolist 
maximises his profit when the buyers represent half of the population. When + 
J) increases above 2 (decreases below 2), the monopolist's optimal price decreases 
(increases) and the corresponding fraction of buyers increases (decreases) . 

Finally, if there are no social effects (J = 0) the monopolist optimal price is a 
solution of the implicit equation: 

1 1 + exp {-(3{pM - h)) 
= Fip,, - h) = • 

The value of l3pM lies between 1 and 1 + exp(/3/i). Increasing (3 lowers the optimal 
price: since the variance of the distribution of willingness to pay gets smaller, the 
only way to keep a sufficient number of buyers is to lower the prices. 

Consider now the case with h < 0, that is, on average the population is not willing 
to buy. Due to the randomness of the individual's reservation price. Hi = H + 9i, 
the surplus may be positive but only for a small fraction of the population. Thus, we 
would expect that the monopolist will maximise his profit by adjusting the price to 
the preferences of this minority. However, this intuitive conclusion is not supported 
by the solution to equations (j25|) when the social influence represented by J is strong 
enough. The optimal monopolist's strategy shifts abruptly from a regime of high 
price and a small fraction of buyers to a regime of low price with a large fraction of 
buyers as (3 J increases. Such a discontinuity may actually be expected for /3J > 4 = 
(3Jb, that is when the demand itself has a discontinuity. But, quite interestingly, 
the transition is also found in the range (3 J a = 27/8 < /3J < 4 == (SJb, that is, 
in a domain of the parameters space (/3J, f5h) where the demand ri{p) is a smooth 
function of the price. 

Such a transition is analogous to what is called a first order phase transition in 
physics 122]: the fraction of buyers jumps at a critical value of the control parameter 
(3Jc{l3h) from a low to a high value. Before the transition, above a value f3J-{f3h) < 
f3JciPh) equations (|^^ already present several solutions. Two of them are local 
maxima of the monopolist's profit function, and one corresponds to a local minimum. 
The global maximum is the solution corresponding to a high price with few buyers 
for pj < [3Jc, and that of low price with many buyers for /3J > 13 Jc- Figure 
presents these results for the particular value f3h = —2, for which it can be shown 
analytically that (3 J- = 4, and (3Jc ~ 4.17 (determined numerically). 

The detailed discussion of the full phase diagram in the plane {f3J,ph}, shown 
on Figure El is presented in the Appendix, and a more general discussion will be 
presented elsewhere [T^ . 
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5 Dynamic features 



In this class of models, the individual threshold of adoption implicitly embodies 
the number of agents each individual considers sufficient to modify his behaviour, as 
underlined in the field of social science j2Ul IH] . We briefly discuss here some dynamic 
aspects, considering a market dynamics with myopic customers: an agent makes its 
decision at time t based on the observation of the behaviour of other agents at time 
t — 1. The adoption by a single agent in the population (a "direct adopter") may 
then lead to a significant change in the whole population through a chain reaction 
of "indirect adopters" lITj. 

Within the McF model, the dynamics for the fraction of adopters in the large N 
limit is then given by 



r]{t) = 1 - F{P -H - J 7]{t- 1)) 



(36) 



and rjit) converges to a solution of the fixed point equation (fT(1|) . As we have seen, 
given J, H and P, two stable and one unstable fixed points appear in (fTU)) for values 
of j3 large enough (small a). The stable solutions correspond to two possible levels 
of r/ at a given price (Figure 0^) . Varying the price smoothly, a transition may 
be observed between these phases. The jump in the number of buyers occurs at 
different price values according to whether the price increases or decreases, leading 
to hysteresis loops. In some cases, the number of customers evolves through a series 
of clustered flips (between = 1 and = 0), called avalanches. For small values 
of (5 (large a), there is a single fixed point for all values of P, and no hysteresis at 

all ( naiiHi). 



1 1,1 l,i 



l,i 1^ 1,E 



a - Hysteresis in the trade-off between prices 
and customers: upstream (biack) 
and downstream (grey) trajectories 




b - Ctironoiogy and sizes of induced effects 
in an avalancfie at tfie pliase 
transition for P = 1.2408 (P„ = 1.25) 



Figure 4: Discontinuous phase transition (full connectivity, synchronous activation regime; source: Phan 
et al. El; parameters: N = 1296, H = I, J = 0.5, /3 = 10). 

The curves in Figure represent the number of customers as a function of the 
price, obtained through a simulation of the whole demand system. The black (grey) 
curve corresponds to the "upstream" (downstream) trajectory, when prices decrease 
(increase) in steps of 10~^, within the interval [0.9,1.6]. We observe a hysteresis 
phenomenon with discontinuous transitions around the theoretical point of symme- 
try, Pn = 1.25. Typically, along the downstream trajectory (with increasing prices, 
grey curve) the externality effect induces a strong resistance of the demand system 
against a decrease in the number of customers. In both cases, large avalanches occur 
at the so-called "first order phase transition". Figure |^ represents the sizes of these 
dramatic induced effects as a function of time (see ^71 for more details). 
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6 Conclusion 



In this paper, we have compared two extreme special cases of discrete choice mod- 
els, the McFadden (McF) and the Thurstone (TP) models, in which the individuals 
bear a local positive social influence on their willingness to pay, and have random 
heterogeneous idiosyncratic preferences. In the McF model the latter remain fixed, 
and give rise to a complex market organisation. For physicists, this model with fixed 
heterogeneity belongs to the class of 'quenched' disorder models; the McF model is 
equivalent to a 'Random Field Ising Model' (RFIM). In the TP model, all the agents 
share a homogeneous part of willingness to pay, but have an additive, time varying, 
random (logistic) idiosyncratic characteristic. In physics, this problem corresponds 
to a case of 'annealed' disorder. In the TP model, the random idiosyncratic com- 
ponent is equivalent to having a stochastic dynamics, because each agent decides 
to buy according to the logit choice function, making this model formally equiva- 
lent to an Ising model at temperature T 7^ in a uniform (non random) external 
field. From the physicist's point of view, the McF and TP models are quite different: 
random field and zero temperature in the McF case, uniform field and non zero tem- 
perature in the TP case. An important result in statistical physics is that quenched 
and annealed disorders can lead to very different behaviours. In this paper we have 
discussed some consequences on the market's behaviour. 

Considering that the seller optimises his own profit, we have exhibited a new 'first 
order phase transition': when the social influence is strong enough, there is a regime 
where, upon increasing the mean willingness to pay, or decreasing the production 
costs, the optimal monopolist's solution jumps from one with a high price and a 
small number of buyers, to one with a low price and a large number of buyers. 

We have only considered fully connected systems: the theoretical analysis of 
systems with finite connectivity is more involved, and requires numerical simulations. 
The simplest configuration is one where each customer has only two neighbours, one 
on each side. The corresponding network, which has the topology of a ring, has been 
analysed numerically by Phan et al. jl7j ) who show that the optimal monopolist's 
price increases both with the degree of the connectivity graph and the range of the 
interactions (in particular, in the case of small worlds). Buyers' clusters of different 
sizes may form, so that it is no longer possible to describe the externality with 
a single parameter, like in the mean field case. Further studies in computational 
economics are required in order to explore such situations. 
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Appendix: Phase Diagram 



In this Appendix we detail the derivation of the phase diagram in the plane {13 J, 
Ph), presented in figure El The phase diagram shows the domain in the parameter 
space where coexist two maxima of the monopohst's profit, one global maximum 
(the optimal solution) and one local maximum. Inside this domain there is a (first 
order) transition line where, as (3 J and/or increases, the optimal solution jumps 
from a solution '-' with a low value = r/_ to a solution '+' with a large value 
r] = rj^, rj being the fraction of buyers. In figure 01 the circles 'o' are points on the 
transition line obtained numerically; all the other curves being obtained analytically 
as explained below. 

In the following without loss of generality we set /3 = 1, which is equivalent to 
say that we measure J and h in units of 1//3. 

To explain the phase diagram in more detail, it is convenient to parametrise 
every quantity /curve as a function of r/. First the {per unit) profit p is given by 

^ JV (37) 



1 — 7] 

(hence the profit A^II = Np{r])rj), and 77 is a fixed point of the equation 

V = G{h,J,r,) (38) 

with 

G{h, J, v) = — . .\j ^ I , (39) 

1 + exp(-/i - 2Jr] + j^) 

We will also make use of an alternative form of (|38j) . (|39j) . that is 

/i = -2j7? + -^ + log(-^) (40) 
1 — r] 1 — r] 

One can also show that the fixed point equation (|38j) is equivalent to 

"'"^'""^O. (41) 

dr] 

and that the condition for having a maximum, H27|l . is equivalent to 

d{r]P'{r])) 



drj 



> 0. (42) 



As we will see, there are two singular points of interest: 
A: Ja = 27/8, /lA = -3/4-log(2); 
B: Jb = 4, hB = -2. 
Let us describe the phase diagram considering that, at fixed J, one increases h 
starting from some low (strongly negative) value. 

If J < 27/8 = J A, the optimal solution changes continuously, the fraction of 
buyers increasing with no discontinuity from a low to a high value as h increases. 
More generally, outside the grey domain in figure |21 there is a unique solution of the 
optimisation of the profit. 
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For J > 27 /8 = Ja, 8iS h increases one will first hit the lower boundary of the 
grey region on the phase diagram, h = h^{J). On this line, a local maximum of 
the profit appears, corresponding to a value r] = r]^ > 1/3. As shown in figure 
Eh , the curve y = G{h, J, rj) intersects y = t] at some small value 77 = and is 
tangent to it at r/ = For h^{J) < h < h+{J) y = G{h,J,r]) has three intersects 
with the diagonal y = r], h = h^{J) being the upper boundary of the grey region 
on the phase diagram. The stability analysis shows that the two extreme intersects 
correspond to maxima of the profit, giving the solutions r] = rj- and r] = r]^. On 
the upper boundary h = h+{J), it is the solution with a small value of rj which 
disappears, with y = G{h, J,rj) becoming tangent to y = 77 for 77 = r/_, see figures 
Etl and \^2. These lower and upper boundaries are obtained by writing that the 
second derivative of the profit with respect to p is zero, giving 

2Jr/(l -r]f = 1 (43) 

Together with (|in|l this gives the curves parametrised by r], 

~ 2r?(l -r/)2 

h = --^ + -^ + log(-^) (44) 
(1 — T])'^ 1 — r] 1 — r] 

the lower curve h-{J) corresponding to the branch r] = r]^ £ [1/3, 1], and the upper 
curve h+{J) corresponding to the branch 77 = r/_ G [0, 1/3]. 

The two curves merge at the singular point A, at which r/+ = 77- = 1/3, Ja = 
27/8, Ha = — 3/4 — log(2). Expanding the above equations (|44j) near r] = 1/3 we find 
that the two curves are cotangent at A, with a slope —2/3. This common tangent 
is thus also tangent to the transition line at A. A straight segment of slope —2/3 
starting from A is plotted on the phase diagram, figure |31 and one can see that this 
is a very good approximation of the transition line for J < 4 = Jb- 

On the lower boundary, for J > 4 = Jb, the local maximum with r] = r]^ appears 
with a negative profit (zero profit at point B where r/_|_ = 1/2). The profit becomes 
positive on the curve starting at point B, on which the profit is zero with ry_|_ > 1/2. 
This curve is obtained by writing p = 0, 77+ > 1/2, that is. 



1 



1 + 



(45) 



and h is obtained as a function of J, by replacing 77 in ()4U() by the above expression 
(j45j) . In this domain of negative profit for the local maximum, the distance to the 
transition line (at a given value of J) is equal to the amount by which the production 
cost per unit of good must be lowered in order to make the solution viable. 

In the domain J > 4 = Jb, the transition line, computed numerically, appears to 
be just above this null profit line. This suggests that an expansion for p small for the 
solution 77+, and for 77 small for the solution ?/_ should provide good approximations. 
The transition is obtained when 11+ = n_ as explained below, and this allows to 
display the curve of figure |31 which turns out to be a very good approximation of 
the transition line for large values of J (or small values of h, typically h < —4.5). 
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Let us first consider the vicinity of the point B at which p+ = 0, r/+ = 1/2, 
and the second derivative of the profit is zero for this '+' local solution. Expanding 
near J = 4, /i just above h-{4) = —2 (p small), one gets the behaviour of the '+' 
solution: 

e = h + 2, < e << 1 
p+ = e + o(e3/2) 

n+ = 1 + 0(6^/2) (46) 

The singular, square-root, behaviour of r/ is specific to point B. For any J > 4, just 
above the null curve r/_|_ increases linearly with e = h — h-{J), but with a similar 
behaviour as for J = 4 for the price and the profit, at lowest order in e: denoting 
by r]^{J) the value r]+{h-{J)), 

< e = /i - /i_(J) << 1 

^+ = ^+^^^ + ^-2Jr;°(l-,0)^ 
P+ = e, 

n+ = r?0(J)e. (47) 

One can see from the expression of r/+ in (|47() how the singularity at point B appears: 
the coefficient of e diverges when condition (|4I'{|) is fulfilled, that is when the solution 
is marginally stable, which is the case at B. 

Similarly one can get the behaviour of the '+' solution near point B at h = —2, 
increasing J from J = 4, as shown in figure |21 




P+ = ^(J-4) + ^(J-4)3/2 

n+ = ^(J-4) + ^(J-4)3/2 (48) 

Coming back to the behaviour at a given value of J, one can get an approximation 
of the '-' solution. The fixed point equation for ij for given values of J and h, is 



77 = H{r]) ^ 1/ 



1 + exp (1 - /i - 2 Jr? + ) 



(49) 



The '-' solution corresponding to a small value of tj can be found by iterating 
ri{k + 1) = H{r]{k)) starting with 77(0) = 0, and r/(/c) is an increasing sequence 
of approximations of ?7_ . The lowest non trivial order is then given by 



77° = iJ(0) = 1/[1 + exp(l - /i)] (50) 
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which is indeed small for h strongly negative. At the next order 



r]i = H{rll) (51) 

Taking rf^ as the small parameter, the expansion of 77! gives 

rt=rt{l + {2J -l)rt) (52) 

and this gives the corresponding approximations for the price and the profit, 

p\. = I -{J- 

= it+Jirtf. (53) 

It is clear from the above equation that the dependency on J is weak since rp_ 
is small, in agreement with the exact behaviour computed numerically, shown on 
figure 12 

Now we consider the neighbourhood of (J, /i_(J)), that \s h = h-[J) + e. De- 
noting by r7o(J) the value of rf^ at h = h-{J), rf_{h) = r]o{J) + er/o(J)(l — r/o(J)). 
Taking this expression for computing the profit of the '-' solution, and writing that 
at the transition the two solutions '+' and '-' give the same profit, one gets the 
following approximation for the value €c{J) of e at the transition (hence the value 
of h at the transition, hc{J) = h-{J) + €c{J))- 

-'^) - ^ <-) 

where r]+{J) is given by equation It is this curve hc{J) = h-{J) + ec(J) which 
is plotted in figure 01 for J > = 4. 
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Figure 5: Functions y = G(/i, J,r7),j/ = r;, price p(t;) and profit n{?7) (H — ). The intersects of j/ = G(/i, J, 77) 
with y = T] give the extrema of the profit; the (possibly local) maxima are those for which dll/dr] > 0. Shown 
here are marginal cases where for one solution dH/dr) = 0, that \s y = G{h, J, rj) is tangent to y = rj (points 
on the lower or upper curves of the phase diagram, figure Isl. 



(a) J = 27/8 h = -3/4 - log(2) 



(h) J = 4: h: 



(cl) Point on the lower curve, 77 < 1/2 (c2) Point on the lower curve, 77 > 1/2 
(dl) Point on the upper curve, J > 4 (d2) Point on the upper curve, J < 4 
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